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SUMMARY 


Two closely related numerical methods which employ operations tables 
have been developed and used in the calculation of the buckling load of 
a monocoque cylinder subjected to pure bending. They are based on the 
assumption of a simplified structure which includes only the most highly 
compressed portion of the cylinder. The first method makes use of a 
lli-row determinant, whereas the second method requires the solution of 
a single 10-row determinant. The buckling loads of three cylinders with 
widely different characteristics were calculated by these methods. 
Reasonable agreement with experiment was obtained. 

A procedure similar to the first method was developed for the cal- 
culation of the buckling load of a cylinder with a cutout. A limited 
experimental check was obtained. 


INTRODUCTION 


The calculation of the buckling loads of reinforced monocoque 
cylinders is a problem of some importance in airplane stress analysis. 
Existing theoretical methods for determining such buckling loads, 
including energy methods, are, in general, lengthy and difficult "to apply. 

A numerical procedure is therefore developed in this report in order to 
simplify the calculations. 

Southwell's relaxation procedure (reference 1) and, in general, methods 
which make use of an operations table (see appendix A) have been success- 
ful in the solution of a variety of stress-distribution problems. It was 
therefore natural that an attempt be made to adapt these methods to 
buckling-load calculations. In reference 2 three closely related methods 
for determining the buckling load from an operations table were established 
and described. A limited experimental check was also obtained. In refer- 
ence 2 the three methods were called the Determinant, Energy, and 
Convergence Methods. In this paper, the first two of these methods, along 
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with slight modifications, are used to calculate the buckling load in 
pure bending of four monocoque cylinders with widely different charac- 
teristics, one of which had a symmetric cutout on the compression side. 

It was found that the buckling loads could be conveniently calculated 
when the actual cylinder was replaced by a simplified structure pre- 
serving the main characteristics of the original cylinder. 

A twofold purpose is thus fulfilled by this investigation. In the 
first place, a method which is fairly short and reasonably accurate is 
developed for the calculation of the buckling load of a monocoque cylinder. 
Secondly, a further experimental check of the methods of reference 2 is 
afforded by a comparison of the theoretical and experimental buckling 
loads for the cylinders considered. 

The authors are indebted to Dr. N. J. Hoff for his advice and help- 
ful criticism, and to Messrs. J. Mele, B. Erickson, and E. B. Beck for 
their part in the experimental phase of the investigation. The work was 
sponsored by and conducted with financial aid from the National Advisory 
Committee for Aeronautics. 


CALCULATION OF BUCKLING LOAD OF CYLINDERS WITHOUT CUTOUT 
Methods of Calculation 


The buckling loads were calculated for three cylinders, the charac- 
teristics of which are given in table I and figure 1. The methods of 
calculation which appeared most convenient are described below. These 
methods yield the load P in the most highly compressed stringer at the 
instant of buckling of the cylinder as a whole. From this load the total 
applied bending moment can be calculated, provided the stress distribution 
is known. The validity and the accuracy of the methods are discussed in 
the next section. Basic theoretical considerations underlying the cal- 
culations may be found in reference 2 and in appendix A. A numerical 
example is given in appendix B. 

Simplified-cylinder solution.- Let the cylinder under consideration 
be replaced by the simplified structure of figure 2. The operations 
table corresponding to this structure is that presented in table II. All 
symbols which appear in this table are defined in appendix C. It should 
be noticed that the operations table is symmetrical about its main diagonal 
The buckling load P has the value which will make the determinant repre- 
sented by table II equal to zero. It may be most conveniently obtained by 
evaluating numerically the determinant for several values of P, plotting 
the determinant values against P, and reading off the load at which the 
determinant is zero. If the load P is lower than the first buckling 


NACA TN 235i| 


3 


load, the determinant will be positive (because it contains an even 
number of rows; see reference 2). If the determinants are evaluated by 
the method of reference 3, that portion of the so-called "auxiliary 
matrix" which corresponds to the first nine rows of table II need be 
considered only once, since it is independent of the load P. 


Solution with assumed displacements .- The above method can be sim- 
plified by assuming the following expressions for the radial displace- 
ments r and the rotations m^ of the most highly compressed stringer 
(stringer 1, fig. 2): 


r = sin^(nx/6L) 


> 

m^ = (dr/dx) = (5>n/6L) sin^(nx/6L) cos (nx/6L) 

J 


( 1 ) 


in which the maximum radial displacement is taken as unity, and L is 
the ring spacing. At rings B, C, and D, respectively (see fig. 2), 
x = L, 2L, and 3L. The determinant is then reduced to that given in 
table III. In the presentation of this table advantage was taken of 
symmetry. The buckling load may be obtained from this determinant by 
the first method given above. Since, however, only the element in the 
lower right-hand corner is a function of the load P, the determinant 
will take the form [K + f(kL)], where K is a constant not dependent 
on P, 


k = 


\|p/(EI ) s t rr 


( 2 ) 


f(kL) is given in table III, and (El) st is the radial bending 

rigidity of the stringer and its effective width of sheet. The value 
of kL at buckling makes the determinant vanish and may be obtained 
from the equation 


f(kL) = -K (3) 

A curve of f(kL) against kL is given in figure 3 and may be used to 
solve this equation in a convenient manner. Consequently the buckling 
load of a cylinder can be obtained from the solution of a single 10-row 
determinant . 

It is useful to note than an upper and lower limit may be found for 
the value of kL at buckling, such that 


1.1*6 < kL < U.ii9 


u 
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The value U.li9 corresponds to the lowest load at which a main-diagonal 
element (in the tenth or eleventh rows) of table II becomes zero. The 
value 1.1*6 is the value at which f(kL) = 0 and is approximate since it 
depends on the assumption of equations (1). 

It may be noticed that both methods require the evaluation of at 
least one determinant. It is suggested that this evaluation be carried 
out by the method of reference 3. The following remarks concerning the 
application of this method in the present problem may be useful: 

(1) The operations table is symmetric about its main diagonal 

(2) The value of the determinant is equal to the product of the 

main-diagonal elements of the auxiliary matrix (defined in 
reference 3) 

(3) The determinant will be equal to zero when the last main- 

diagonal element of the auxiliary matrix vanishes (see 
appendix A) 

Discussion of Methods 

The methods outlined in the preceding section are based on the 
simplified structure of figure 2. The following considerations underlie 
the choice of this structure and of the methods of calculation: 

(1) The most highly compressed stringer was considered of paramount 
importance at buckling, so that it was thought permissible to neglect 

all other stringers in these approximate calculations. This is equivalent 
to considering the most highly compressed stringer as a column elastically 
supported by the rings and sheet. The elasticity of the supports is 
represented by the ring and sheet influence coefficients in the operations 
table (appendix B). 

(2) Points on the tension side of the cylinder will undergo only 
negligible displacements and hence may be considered fixed. The rings 
are therefore assumed to continue up to a point, near the tension side, 

90° away from the most highly compressed stringer, and to be rigidly 
fixed there (fig. 2). 

(3) It would seem natural to continue the sheet up to the same point 
as the rings. Because all stringers except the most highly compressed 
one have been neglected, this would imply a single panel of sheet in 
each bay, extending over 90°. The operations table, however, is set up 
considering each panel with its edge reinforcements as a unit in which 
only the corner points have independent freedom of motion (see, e.g., 
reference U). Therefore the action of the 90° sheet panel would be 
determined by the displacements of its corners, with no possibility of 
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intermediate adjustment. Consequently the rigidity of the panel would 
be greatly exaggerated. The decrease in the effective shear modulus of 
the buckled sheet (reference 5) because of the larger angle subtended 
would not provide a sufficient reduction in the shear rigidity. The 
sheet panel was therefore taken to be smaller, the natural stopping point 
being the position of the stringer next to the most highly compressed 
stringer in the actual cylinder. Thus only the sheet which provides 
additional stiffness to the most highly compressed stringer is considered. 
This appears consistent with the assumption that all other stringers may 
be neglected. A point with independent freedom of motion was therefore 
considered in each ring at the intersection with the free edge of the 
sheet. It may be remarked that, if the rings were to be terminated 
there, the consequent reduction in the influence coefficients would in 
general be negligibly small. 

(U) The length of the cylinder was considered constant and equal to 
six times the ring spacing. For the three cylinders investigated, PIBAL 
cylinder 10 and GALCIT cylinders 25 and 65 (fig. l), this corresponds to 
1.5, 1.5, and 1.2 times the respective diameters. For the fuselage of a 
large modern transport this length would be approximately equal to the 
diameter. The following table may be set up on the basis of experimental 
results presented in the references given: 


Loading 

Limiting 

Increase in buckling load 


value of 

L'/D 

at lower value of L'/D 
(percent) 

Reference 

Compression 

1.5 

6 at L'/D = 1.0 

6 

Pure bending 

2.0 

12 at L'/D = 1.2 

7 


The buckling load is practically independent of the length if the length- 
to-diameter ratio L'/D is equal to or larger than the limiting value 
given. Examples of the increase in the buckling loads for cylinders 
shorter than the limiting length are given in the third column of the 
above table. The length assumed in the calculation will be in general 
somewhat shorter than the limiting length; the error caused by this may 
be estimated with the aid of the above table to be at most 10 or 15 per- 
cent of the buckling load of a cylinder longer than the limiting length. 
The effect of the length was investigated in some detail with test 
cylinder 25 of the GALCIT series (reference 8). The buckling load for 
this cylinder was calculated considering different numbers of bays and 
the results are shown in figure i;. It may be seen that the calculated 
buckling loads approach some constant value in what appears to be an 
asymptotic variation and that the difference in the buckling loads 
obtained considering six or eight bays is small. It was concluded that 
the small improvement in accuracy given by a longer structure did not 
warrant the increased amount of work required to obtain it. 
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(5) The buckling load was calculated with the aid of the simplified 
structure for test cylinders 25 and 65 of the GALCIT series (reference 8) 
and for cylinder 10 of the PIBAL series (reference 9). Those specimens 
were chosen because of their widely different characteristics (see fig. 1 
and table I). Comparisons of the results of the present analysis with 
those of experiment are presented in table IV and in figures 1;, 5, and 6. 
The calculated buckling loads may be seen to be consistently higher than 
the corresponding experimental values. The percentage errors obtained 
are not considered excessive, however, upon comparison with the results 
obtained earlier at PIBAL by means of strain-energy methods. One of those 
solutions (reference 10) gave better results than the present investi- 
gation, but required a prohibitive amount of work. 

(6) Approximate deflected shapes at buckling obtained with the aid 

of the simplified structure of figure 2 are given in table V for the three 
cylinders investigated. The same table also gives results of measurements 
made on some actual test specimens after buckling (reference 6). It may 
be noticed that fair agreement has been obtained between measured and 
calculated values, so that an additional indirect experimental check has 
been provided on the reasonableness of the simplified structure chosen. 

It should be kept in mind that the measurements were taken after the 
cylinders had buckled, and therefore may differ from the actual displace- 
ments at the instant of buckling. 

(7) Table V also shows that the radial displacements r of the most 
highly compressed stringer at buckling are closely represented by 

r = sin n (nx/6L) (1*) 

where n = 1*, 5, or 6. The rotations m^ may be closely approximated by 

m-j. = (nn/6L) sin n “^(nx/6L) cos (nx/6L) (5) 

The values obtained with n = 5 represent a reasonable average of all 
experimental and calculated deflections, and therefore this value of n 
was chosen for the solution with assumed displacements which was described 
previously. As a check, the buckling load of GALCIT cylinder 65 was cal- 
culated by that method and was found to be 173Q pounds, the buckling load 
calculated from the simplified structure without the assumption of dis- 
placements was 1670 pounds (fig. 5), so that the error introduced by this 
assumption is only 3.6 percent of the latter value. Cylinder 65 was 
chosen since it is the specimen for which the agreement between assumed 
and actual displacements is the poorest (table V). It should be remembered 
in this connection that it was shown elsewhere (reference 2) that the 
methods of calculation used in this report are not too sensitive to errors 
in the assumed deflected shape. 


NACA TN 235U 


7 


(8) The operations tables for the simplified structure (tables II 
and III) may be put in nondimensional form by the following process. 
Divide the tenth and eleventh rows and columns by L, and all terms in 
rows and columns (10) through (li|) by the quantity G e fftd/L. The beam- 

column terms appearing in the lower right-hand corner may then be written 
as (El) stri> ^L^ times some function of kL. The ring influence coef- 
ficients rr, rm, rt, and so forth are equal to (El) r /d? multiplied 
by some function of r/d (reference 11). If all rows and columns are 
divided through by (El) r /d^ it will be noticed that the operations 
table will be a function of the four nondimensional parameters 


r U ( EI )str r 

< EI )r 


r 


Gefft* 

(EDji 

r/d 


(6) 


kL 

where r is the cylinder radius, (El) r the bending rigidity of a ring 
in its own plane, d the circumferential stringer spacing, G e ff the 

effective shear modulus for the sheet, t the sheet thickness, and the 
other symbols have been previously defined. The effects of shearing and 
extensional deformations of the rings, respectively, are represented by 
the two additional parameters: 


£ = A r */A r 
7 = A r d 2 / l r 


where I r is the moment of inertia of the ring cross section, A r is 
the area of the ring, and A r * is the effective shear area of the ring 
cross section. Reference 11 shows, however, that the effect of these two 
parameters is in general negligible. 

It has been shown in reference 12 that the buckling load of a mono- 
coque cylinder depends on the parameter A. Two additional parameters 
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were established in that reference by physical reasoning to be r/d and 
e/€ cr , where e is the strain in the most highly compressed stringer 

at failure and e is the buckling strain of a sheet panel. An experi- 
c r 

mental verification of the fact that these parameters approximately 
control the buckling phenomenon in monocoque cylinders is given in refer- - 
ence 13. It may be seen that two of the parameters found in the present 
development are the same as those found in reference 12, while the param- 
eter T includes the quantity eA cr , since the shearing rigidity G eff 

was found in reference 5 to be closely approximated by 


^ = N + (1 - N)e N €/€cr (8) 

G 0 

where 

N = 0.0275 [(2nr/d) + l] (8a) 

and Gq is the shear modulus of the sheet material. 

Curves of kL, which represents the buckling load, plotted against 
the parameters A and r are shown in figure 7 for all the cylinders of 
reference 8 with r/d = 6.32. An insufficient number of cylinders is 
available so that the position of these curves is not definitely determined. 
It may be stated, however, that the results presented do not contradict 
the validity of the four parameters established. 


CALCULATION OF BUCKLING LOAD OF A CYLINDER WITH A CUTOUT 
Experimental Investigation 


The methods developed previously were extended to include cylinders 
with cutouts. A cylinder with a cutout (PIBAL cylinder 82, fig. 8), which 
consisted of a thin circular shell reinforced by six stringers and four 
evenly spaced rings was therefore constructed and tested. The cutout 
extended circumferentially for 90° on the compression side of the cylinder. 
Pure bending moments were transmitted to the ends of the cylinder through 
heavy rings which could be assumed rigid. The test rig was the same as 
that used in the cylinder tests of reference 13. This cylinder buckled 
when the load in the most highly compressed stringer (stringer 2 in fig. 8) 
was 5U00 pounds. This load corresponded to a total applied moment of 
158,000 inch-pounds. Photographs of the buckling cylinder are shown in 
figures 9 and 10. 
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Method of Calculation 

The results of the theoretical investigation indicated that the 
simplified methods evolved for complete cylinders are not satisfactory 
for cylinders with cutouts. For such cylinders, the distortion under 
load depends primarily on the geometry of the cutout, so that only dis- 
placements in the vicinity of the cutout require consideration in the 
operations table . 

With consideration of the symmetry of the cylinder, an operations 
table including only the displacements r B1 and r B 2 and the rotation 

mj.g 2 (see fig. 8) was set up. The resulting buckling load was, for all 

practical purposes, the same as that obtained through the use of an 
operations table which permitted all possible generalized displacements 
at all the joints of the cylinder. Hence the assumption that only the 
joints in the vicinity of the cutout need be considered in the operations 
table was justified for the cylinder with a cutout. The actual experi- 
mental deflected shape of PIBAL cylinder 82 (figs. 9 and 10) shows that 
the major distortions took place in the vicinity of the cutout, and that 
all other joints may be assumed to have had zero displacements. 

Cylinders encountered in practice, however, will be of a more compli- 
cated construction than PIBAL cylinder 82, and hence the simplified 
operations table described above may not be sufficiently complete. 
Depending on the size of the cutout, it is suggested that the operations 
table be expanded so as to include all the joints surrounding the par- 
ticular cutout. 

The buckling load obtained for PIBAL cylinder 82 considering only 
three generalized displacements was 81*00 pounds. The discrepancy between 
the theoretical and experimental buckling loads was attributed, mainly, 
to the inaccuracy of the value of the effective shear modulus Ggff used 

in the calculations. This value was 0.71G 0 and was taken from equa- 
tion (8). This equation is based on bests on panels buckled because of 
compression. The sheet panels in the present cylinder, however, are 
under the action of combined compression and shear. No values for the 
effective shear modulus of curved panels under such a loading could be 
found in the literature, but, according to data obtained from flat panels 
(reference ll*), it appears that the correct value of G e ff should be 

considerably lower. Furthermore, as is shown in the next section, there 
is reason to reduce the shear modulus even further. 

The calculations for PIBAL cylinder 82 were therefore repeated with 
an assumed value of G e ff = O.IGq. The resulting buckling load was 

5900 pounds, which may be seen to be in good agreement with experiment. 
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Reduced Effective Shear Modulus 

If a panel of sheet is not in a buckled state, the relation between 
shear stress t and shear strain 7 is simply Hooke's law: 

'r = (?) 

If the panel is in a buckled state, a relation analogous to equation (9) 
will still hold between the average shear stress x„ ir and the average 

shear strain 7 av > provided that an effective shear modulus G e ff is 
used in place of' Gq. In other words, 

T av ~ G eff^av ( 10 ) 

and G eff " G 0 if the panel is not buckled. The value of G eff will 

represent the complex state of stress of the buckled panel, and will 
presumably vary with panel dimensions and type of loading. 


The value of G e ff is the proportionality factor between the average 

shear stress and the average shear strain. In problems of instability, 
however, it is desired to know the relation between a small increase in 


stress d 


( T av) 


and a small increase in strain 




This relation 


will again have the same form as equation ( 9 ), if only a reduced effective 
shear modulus G e ff rec } is use d i n place of Gq. In other words, 


d ( T av) G eff red d ( y av) 


( 11 ) 


Thus this new modulus represents the resistance the panel will offer 
against distortions additional to those represented by y a y According 

to the previous discussion, this new modulus will also depend upon the 
dimensions of the panel and upon the amount of shearing and compressive 
loads present. 


If equation (10) is written in differential form as 


^ T av) 


7 


d ( G eff ) + 


av 


G eff 


d hav) 

comparison with equation ( 11 ) indicates that 

d ( G eff ) 
d ( y av) 


d ( 7 av) 


G eff r ed = G eff + y av 


( 12 ) 


(13) 
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The following remarks may be made about the reduced effective shear 
modulus G e ff 

(1) Geff re d = Go when the panel is not in a buckled state 

(2) G effred = G eff when the average shearing strain in the panel 
is zero immediately before buckling 

(3) G e ff red = G e pf when there is no change of shearing strain 
during buckling of the structure under consideration 

(U) G e pp red < G e ff in all other cases, since in general the modulus 
Geff decreases with increasing shear strain y av , so that the second 
term in the right-hand side of equation (13) is negative 

The latter case applies to the cylinder with a cutout. The low 
value assumed for the shear modulus in the calculations is therefore 
plausible . 


Polytechnic Institute of Brooklyn 

Brooklyn, N. Y., August 31, 19U8 
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APPENDIX A 
BASIC THEORY 


The procedures developed in this report for the calculation of the 
buckling load of a monocoque cylinder are based on methods developed in 
reference 2 which make use of an operations table similar to that used 
in Southwell's method of systematic relaxation. These methods are out- 
lined here, more rigorous proofs being given in reference 2. Rigorous 
proofs are only given here for some modifications of these methods which 
were not discussed in that reference. 

Consider several points in the structure in question distributed so 
as to cover the entire structure. Let these points be numbered consecu- 
tively from 1 to n. The generalized force exerted on joint i by a gener- 
alized displacement xj at joint j (all joints but j being considered 

temporarily rigidly fixed) may be denoted by a-j^Xj. The quantity a-y 
is called an influence coefficient. If Fj^ is the generalized external 
force acting at joint i, the equilibrium condition for the ith joint is 

n 

F i + XZ a ij X j = ° (A1) 

j= 1 

provided that the principle of superposition is valid. If equation (Al) 
is written for every joint in the structure, a set of linear simultaneous 
equations will result with generalized displacements as unknowns. The 
array, or matrix, of the coefficients of this set of equations is called 
the operations table and may be written as 


a ll 

a 12 * 

. . a H . 

. . a lj . 

• • a ln 

a 21 

• • 

a 22 • 

. . a 2i . 

. .. a 2 j . 

. . a 2n 

• • • • 

• • 

a il 

a i2 • 

. . a ii . 

. . a ij . 

• • • • 

• • a in 

a jl 

• • • 

a j2 • 

. . aji . 

* * a jj • 

• * a jn 

• • 

a nl 

a n2 • 

. . a ni . 

. . a n j .' 

• • a nn 


As a consequence of Maxwell 1 s reciprocal theorem aj_j = a ji 
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Equation (Al) represents the equilibrium conditions for the given 
structure in terms of displacements. In general, the determinant A n 

will not be equal to zero; then, only one set of displacements may be 
found which will satisfy the equilibrium conditions. If the determinant 
A n vanishes, however, more than one such set of displacements will exist. 
This is physically possible only at neutral equilibrium, or, which is 
the same, at a buckling load. This leads to what was called in refer- 
ence 2 the Determinant Method, the basis of which is the fact that the 
lowest load at which the determinant A n vanishes is the lowest buckling 
load. 


A proof will now be given of the fact that in general at the lowest 
buckling load the last main-diagonal element of the auxiliary matrix of 
the method for evaluating determinants given in reference 3 is equal to 
zero. Let the symbol A± stand for the determinant 


A i 


a ll a 12 • * • a li 
a 21 a 22 • • • a 2i 


a il a i2 • • • a ii 


(A3) 


The value of this determinant is equal to the product of the first i main- 
diagonal terms of the auxiliary matrix. If auxiliary-matrix elements are 
denoted by the symbol aj_ j , then 

i 

A i = IT a jj = a ll a 22 • • • a kk • • • a ii (A b) 

j=l 


Let A^ be the first of these determinants to vanish; then theorem 2 
of reference 2 gives 


A k ~ A k+1 - • • • ~ A n - 0 (A5) 

where A n is the determinant given in equation (A2). Two cases may 
then be considered: 

Case 1; k = n .- In the case where k = n, A n is the only one of 
these determinants which vanishes. By equation (AJU) the only factor 
which is contained in A n and in no other A^_ determinant is o^, 

which, as was to be proved, must therefore vanish. 

Case 2; k < n .- In the case where k < n, the method of reference 3 
fails to give any terms beyond a^, which of course is zero. Here a m 
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u* 


is obviously not the first term to vanish, but in this case several of 
the higher buckling loads are identical with the first. This case is 
expected to occur rather infrequently. 

The Energy Method of reference 2 is based on the condition that the 
second variation of the total potential energy must vanish at buckling. 
This condition may be written as 


This equation is satisfied at buckling by the buckling displacements. 

In the Energy Method some of these displacements, say Xp, Xp+j, • • • , 
x n , are guessed; then the others are obtained from the conditions 


The matrix of the coefficients of these simultaneous equations, including 
the constant terms, is the reduced matrix A', where: 



|£L =0 k = 1, 2 


. ., p - 1 


(A7) 


a ll a 12 * * * a l,p-l a 'l,p 

a 21 a 22 • * • a 2 , p-1 a '2,p 


A' 


(A8) 



in which 


n 



J=P 


and 


n n 



i=P j=P 


It will now be proved that the Determinant Method may be applied to the 
operations table of equation (A8). 
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It is well-known (see reference 15) that any quadratic form Q may be 
put in the form 


Q = 


t 

n =1 


bi(Li) 


2 


(A9) 


where the bj_ quantities are constants, and 

n 

lq = c j x j (AlO) 

j=l 

where the cj quantities are constants. It is assumed that there are 
n linearly independent quantities L i . This assumption entails no loss 

of generality since in the case in which it is not true some of the 
constants b-^ will be zero. 


By means of equation (A 9 ) and table 2 of reference 2 the following 
table may be set up: 


Sign of bj_ 

Sign of Q 

Classification of Q 

Type of 
equilibrium 

All, less 
than zero 

Q < 0 always 
(Q = 0 if xj = 0) 

Negative definite 
nonsingular 

Stable 

Varying ; none , 
zero 

Q may be positive, 
negative, or zero 

Indefinite 

nonsingular 

Unstable 

All, greater 
than zero 

Q > 0 always 
(Q = 0 if xj=0) 

Positive definite 
nonsingular 

Unstable 

Some, zero; all 
others, negative 

Q = 0 or 
Q < 0 

Negative definite 
singular 

Neutral 

Some, zero; 
others, varying 

Q may be positive, 
negative, or zero 

Indefinite singular 

Unstable 

Some, zero; all 
others, positive 

Q = 0 or 
Q > 0 

Positive definite 
singular 

Unstable 
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Two conditions for neutral equilibrium have been thus set up: 

(1) The vanishing of the determinant A n of the quadratic form Q 

(2) The fact that some of the b^'s equal zero, while all others 
are negative 

As both conditions are necessary and sufficient, they are equivalent and 
may be used interchangeably. 

When some of the displacements are guessed as previously explained, 
the quadratic form Q becomes the reduced quadratic form Q f , where 

p 

Q' = H bi'(Li') 2 

i=l v 

where the bj_' quantities are constants and 

^i 1 = c p' + y - c j ' x j 
3=1 J 

The statements of the above . table may now be applied to the quadratic 
form Q', since Q' is the value the quadratic form Q will take on 
when the displacements Xp, . . ., x n . are assumed. Neutral equilibrium 
will then exist when some of the b^'s are zero and all others are neg- 
ative. But this condition is equivalent to the condition that the determi- 
nant corresponding to the quadratic form Q' must vanish. This determinant 
is obtained by multiplying out the right-hand side of equation (All) and 
expressing the x-esult in the form 


P-1 P-1 

P-1 

P-1 


«• 

a ij x i x j a 'ip x i 

+ ^ a 'pi x i + a 'pp 

(A13) 

i=l 3=1 

i=l 

3=1 



The determinant corresponding to Q' will then be seen to be identical 
with A' of equation (A8). 

It therefore follows that the vanishing of the determinant A’ of 
equation (A8) corresponds to neutral equilibrium. If the displacements 
Xp, • • • , x n were not chosen exactly equal to the displacements of 

the structure at buckling, an approximate value of the buckling load will 
be obtained rather than an exact one. It was proved in reference 2 that 
this approximate load will be higher than the actual one. 


(All) 


(A12) 


3 
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APPENDIX B 


GENERAL FORMULAS AND NUMERICAL EXAMPLE 


In this appendix is presented the procedure for the determination 
of the influence coefficients required in setting up the operations 
tables II and III. Since many of the formulas used in the analysis are 
scattered throughout the literature, some of these are given here, 
together with appropriate reference. When a formula is not listed, refer- 
ence to its source is given. A numerical example is also given illus- 
trating both methods suggested for the calculation of the buckling loads 
of cylinders without cutout. 


The operations tables (tables II and III) contain three types of 
influence coefficients, which represent the effects of the rings, the 
sheet covering, and the stringers. The methods by which each type is 
determined are outlined below. 

Ring influence coefficients .- The ring influence coefficients are 
characterized by the symbol as, for example, rr^, rt^, or tnp. 

These coefficients may be determined from reference 11, in which they 
may be seen to depend on the three parameters p, y, and %, where p 
is the central angle of the ring segment, and y and £ are defined in 
equations (7). Any one of the three following ways may be used for the 
calculation of these coefficients: 

(1) General formulas are given in equations (20), (27), (28), and 
(29) of reference 11. As an example, the formula for rtpj is repeated 
here with a slight change in notation: 


Influence Coefficients 



(Bl) 


where 


-A tr pU(Ei) ftrl + \\ f tr2 


I 
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and 


,2 P 7 


= fil + (1/y) A2 + '2/ 0 (f £3)j + (l/?) 2 < (1/^) + 3 1/ £) - 7j 2 ^aU 


in which the quantities f are functions of the central angle p and 
are given in reference 11 by equations ( 2h ) and, for specific values 
of p, by table II and figures 3 to 13. 

(2) Influence coefficients obtained from the formulas mentioned 
under item (l) are tabulated in tables III and IV of reference 11 for 
specific values of the parameters p, y, and I. 

(3) From the tabulated values mentioned in item (2) above, curves 
were plotted which are presented in figures llj. to 85 of reference 11. 
It should be noted that in general these curves are accurate only for 
values of p = 15 ° • 

Sheet influence coefficients.- The shear in the sheet covering is 
represented by the influence coefficients containing the quantities a 
and Aj, as, for example, 2a r 2 Aj or 2a t a n Aj dj. The quantities a 

may be determined in good approximation from the formulas 




a r = O.ip 

04 . = - 0 . 5(1 - 0.01666 • • • p 2 ) 

a n = -0.008333 •••?(!+ O.OlU 286 p 2 ) 




(B2) 


taken from page 27 of reference 11. It should be noted that in this 
reference a r , a t , and a n are denoted as r q ^(Lq), tq^(Lq), and 

n q/(l 2 q), respectively. 

The quantity Aj is given by 


A x = 


I ~ G eff 


tdx 


(B3) 


where the effective shear modulus Gg-ff is given by equation (8) of the 
present report. 
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The buckling strain e cr of a sheet panel appearing in this equa- 
tion was obtained by means of Redshaw's formula: 


where 


e cr " ( e flat/ 2 ) + \J( e flat/ 2 ) 2 + ( e curved) 2 


(B h) 


c flat 


k'n 2 /t\ 2 
12(l - V 2 ) 


(B5) 


as shown in reference 16, in which k' is the end-fixity coefficient, 
V is Poisson's ratio, and 


e curved- 


0.6(t/r) — ~ hi x 10 7 (r/t) 2 
1 + O.OOh (E/F cy ) 


(B6) 


as shown in reference 17, in which F cy is the yield stress of the 
material . 


The effect of normal stresses in the sheet covering is taken into 
account by an effective width of sheet as discussed in the next section. 

Stringer beam-column influence coefficients .- The beam-column 
effects in the stringers are represented by the influence coefficients 
containing the load P, as, for example, (p/Dq)(l - c) or Pks/Dq. 

A complete list of beam-column influence coefficients is presented in 
tables VI, VII, and VIII in which the sign convention as well as defi- 
nitions of symbols are given. All these coefficients are functions of 
the quantity k defined in equation (2). The following formula which 
was derived in reference 18 on the basis of work contained in reference 19 
is suggested for the calculation of the effective width 2w: 


2w = (1/ e ) (d/r) Jo.3t + 1.535 [(t/d) (re 



(B7) 


If the load P causes a stringer stress which is higher than the 
proportional limit of the material, the modulus of elasticity £ which 
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appears in equation (2) must be reduced in an appropriate manner. In 
this investigation Von Karmen's formula was used: 


E red = h EE-b/(vlE + vjE t ) 2 (B8) 

where is the tangent modulus. 


Numerical Example 

As an example of the application of the methods suggested for the 
calculation of the buckling load of cylinders without cutout the buckling 
load is determined here for GALCIT cylinder 65 of reference 8. The 
characteristics of this cylinder are given in figure 1. Table I is 
readily set up with the aid of these characteristics and the equations 
listed in the previous section. The influence coefficients required 
for tables II or III can then be calculated. 

Ring influence coefficients.- The values of the ring influence 
coefficients corresponding to Pj, y j, and lj (see table I) happen 

to appear in tables III of reference 11. However, since these tables 
do not contain coefficients corresponding to Ijj, it is necessary to 

use values interpolated from the appropriate curves of reference 11. 

The following values were obtained for the ring influence coefficients: 


Ring segment I 


■ 6 - 393 


■ - 20 - 38 


’ 109 -° 


S x® ' ?8.U3 


Ring segment II 


Sills Tr ■ 8 - 65 


rn. 


( d II 2 
Mill EI r 


= -33.0 


tni 


it 


ii 


lM Il\ Elj 


= 33.0 


- 170 
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‘AE ’ - 656 ' 6 


^iifer '-- 200 


ttM i iir ■ U999 


tt M . 


Wn 3 

II\ EI„ 


280 


rn,- 


'I El, 


= - 8.^31 


rr F, El? * ^' 86 


Msfc 1 " 659 - 7 


Sheet influence coefficients .- The quantities a of equation (B2) 
have the following values: 


Bay II 

No sheet in this bay 


Bay I 

a r = 0.0261 

0 ^ = -0.1+99 

an = -0.00218 

Stringer beam-column influence coefficients .- For an assumed load 

P = 1550 pounds a typical set of calculated results is given below (see 
tables I and VI and equations (2) and (B8)). 


p 

(lb) 

a 

(psi) 

Ered 

(psi) 

Ered ^str r 

k 

kL 

sin kL 

cos kL 

1560 

37,300 

7.056 x 10 6 

1+000 

0.6225 

2.1+90 

0.6061+6 

- 0.79511 
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The beam-column influence coefficients required for table II are 

then 


^ = 281.31 P ( k - k~ s . ) = 225k. 6 

kDq 

^ = 1337.6 (s - kLc) = 6191.7 

No beam-column influence coefficients need be calculated if table III 
and figure 3 are used. 

Calculation of buckling load by table II.- The state of stability 
of the structure at a load of P = 1550 pounds may now be investigated 
by introducing all the above influence coefficients into table II and 
evaluating the corresponding determinant. The results corresponding to 
the assumed load of 1550 pounds, as well as those for the loads of 1650, 
1680, and 1715 pounds, are presented in figure 5 . In this figure are 
plotted both the values of the determinant and the values of the last 
main-diagonal term a nn of the auxiliary matrix of reference 3- The 
intercept of the curves in this figure may be read off and corresponds 
to the buckling load. It should be noted that only the stringer beam- 
column influence coefficients vary as the assumed load is changed. 

Calculation of buckling load by table III.- The value of the 

quantity K of equation (3) was found by the method of reference 3 to 
be -5.1*72. From figure 3, this value may be seen to correspond to 
kL = 3.1*, from which P cr = 1730 pounds. Care must be taken that a 

reduced modulus (equation (B8)) be used if the stress at buckling is 
above the proportional limit. 
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^1’ ^2’ 

Ai, • * • A n 

^ e ^str 

A r 

A r * 

A' 

A, B, C, D 
E 

(El) r 

( £I ) S tr r 

E red 

E t 

F cy 

Fi 

G eff 

Ge f fped 
G 0 

Ir 


APPENDIX C 
SYMBOLS 

determinants or matrices 

effective cross-sectional area of stringer 
effective cross-sectional area of ring 
effective shear area of ring cross section 

reduced matrix 
rings 

Young 1 s modulus 

bending rigidity of a ring in its own plane 

radial bending rigidity of stringer and its 
effective width - of sheet 

reduced modulus 
tangent modulus 
yield-point stress 
generalized external force 
effective shear modulus 
reduced effective shear modulus 

shear modulus of sheet material under no com- 
pressive load 

moment of inertia of ring cross section plus 
effective width 
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2li 

•^strj, 

K 

L 

L i 

L'/D 

Mt 

N 

P 

Per 


Q 

Q' 

R 

T 



b±, V 

°j’ V 

d 

f 

f(kL) 

i, j 


moment of inertia of stringer plus effective width of 
sheet for radial bending 

constant 

ring spacing 

linear function 

ratio of total cylinder length to cylinder diameter 

moment causing bending of stringer (vector pointing 
in tangential direction) 

moment causing bending of ring in its plane; also 
0.0275 [(2nr/d) + l] 

axial stringer load 

load in the most highly compressed stringer at the instant 
of buckling 

quadratic form 

reduced quadratic form 

radial force 

tangential force 

element of operations table 

element of last row or column of reduced matrix 

constants 

constants 

circumferential stringer spacing 
function of |3 
function of kL 
indices 


k =\|P/( EI ) s tr r 
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v « 

k' 

end-fixity coefficient 

H 

rotation causing bending of stringer (vector pointing 
in tangential direction) 

n 

number of stringers; number of generalized displacements; 
rotation causing bending of ring in its plane 

p 

index 

r 

radius; radial displacement 

/"-S 

rr, rm, rt, 
and so forth 

ring influence coefficients 

t 

sheet thickness; tangential displacement 

2w 

effective width of sheet 

X 

longitudinal axis of stringer 

x j 

generalized displacement 

1 

r 

parameter 

K 

parameter 

A i 

/Ggfftd-A 
parameter ( 1 

a ij 

auxiliary-matrix element 

°nn 

last main-diagonal element of auxiliary matrix 

a t’ a n 

functions of j3 required for sheet influence 
coefficients 

P 

central angle of a ring segment (d/r) 

7 

_ 

parameter; shear strain 

e 

strain in most highly compressed stringer at failure 

€ cr 

buckling strain of a sheet panel 

. ^curved 

buckling strain of nonreinforced circular cylinder 
under uniform axial compression 
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e flat 

v 

I 

CT 

T 

Subscripts : 


buckling strain of flat panel under uniform compression 

Poisson's ratio 
parameter 

compressive stringer stress 
shear stress 


A, B, 

C, D 

rings 

F 


fixed 

M 


movable 

av 


average 

exp 


experiment 

1, 2, 

3 

stringers 


I, II 


regions 


NACA TN 2351i 


27 


REFERENCES 


1. Southwell, R. V.: Relaxation Methods in Engineering Science. A 

Treatise on Approximate Computations. The Clarendon Press (Oxford) 

19l*0. 

2. Boley, Bruno A.: Numerical Methods for the Calculation of Elastic 

Instability. Jour. Aero. Sci., vol. ll*, no. 6, June 19U7, 

PP- 337-31*8 . 

3. Crout, Prescott D.: A Short Method for Evaluating Determinants and 

Solving Systems of Linear Equations with Real or Complex Coefficients. 
Trans. AIEE, vol. 60, 19l*l, pp. 1235-121:0. 

1*. Hoff, N. J., Levy, Robert S., and Kempner, Joseph: Numerical Pro- 

cedures for the Calculation of the Stresses in Monocoques. I - Dif- 
fusion of Tensile Stringer Loads in Reinforced Panels „ NACA 
TN 931*, 19l*U. 

£. Hoff, N. J., and Boley, Bruno A.: The Shear Rigidity of Curved Panels 

under Compression. NACA TN 1090, 192*6 . 

6. GALCIT: Some Investigations of the General Instability of Stiffened 

Metal Cylinders. IV - Continuation of Tests of Sheet-Covered 
Specimens and Studies of the Buckling Phenomena of Unstiffened 
Circular Cylinders. NACA TN 908, 19l*3. 

7. GALCIT: Some Investigations of the General Instability of Stiffened 

Metal Cylinders. VI - Stiffened Metal Cylinders Subjected to 
Combined Bending and Transverse Shear. NACA TN 910, 19 ^ 3 . 

8. GALCIT: Some Investigations of the General Instability of Stiffened 

Metal Cylinders. V - Stiffened Metal Cylinders Subjected to Pure 
Bending. NACA TN 909, 191*3. 

9. Hoff, N. J., Fuchs, S. J., and Cirillo, Adam J. : The Inward Bulge 

Type Buckling of Monocoque Cylinders. II - Experimental Investi- 
gation of the Buckling in Combined Bending and Compression. NACA 

TN 939, 191*1*. 

10. Hoff, N. J., Klein, Bertram, and Boley, Bruno A.: The Inward Bulge 

Type Buckling of Monocoque Cylinders. V - Revised Strain Energy 
Theory Which Assumes a More General Deflected Shape. NACA TN 1505 
191 * 8 . 


28 


NACA TN 235U 


11. Hoff, N. J., Klein, Bertram, and Libby, Paul A.: Numerical Pro- 

cedures for the Calculation of the Stresses in Monocoques. IV - 
Influence Coefficients of Curved Bars for Distortions in Their 
Own Plane. NACA TN 999, 19U6. 

12. Hoff, N. J. : General Instability of Monocoque Cylinders. Jour. 

Aero. Sci., vol. 10, no. h, April 19h3, pp. 105-llU, 130. 

13. Hoff, N. J., Boley, Bruno A., and Nardo, S. V.: The Inward Bulge 

Type Buckling of Monocoque Cylinders. IV - Experimental Investi- 
gation of Cylinders Subjected to Pure Bending. NACA TN 1U99, 

191+8 . 

111. Kromm, A., and Marguerre, K.: Behavior of a Plate Strip under 

Shear and Compressive Stresses beyond the Buckling Limit. NACA 
TM 870, 1938. 

15. Bocher, Maxime: Introduction to Higher Algebra. The Macmillan Co., 

1931. 

16. Timoshenko, S.: Theory of Plates and Shells. First Ed., McGraw- 

Hill Book Co., Inc., 19l|0. 

17. Donnell, L. H.: A New Theory for the Buckling of Thin Cylinders 

under Axial Compression and Bending. Trans. A.S.M.E., vol. 56, 
no. 11, Nov. 193b, PP* 795-806. 

18. Hoff, N. J., and Klein, Bertram: The Inward Bulge Type Buckling 

of Monocoque Cylinders. I - Calculation of the Effect upon the 
Buckling Stress of a Compressive Force, a Nonlinear Direct Stress 
Distribution, and a Shear Force. NACA TN 938, 19U1|. 

19. Ebner, H.: The Strength of Shell Bodies - Theory and Practice. 

NACA TM 838, 1937* 

20. Templin, R. L., Hartmann, E. C., and Paul, D. A.: Typical Tensile 

and Compressive Stress-Strain Curves for Aluminum Alloy 2i|S-T, 
Alclad 2J|S-T, 2i|S-RT, and Alclad 2i|S-RT Products. Tech. Paper 
No. 6, Aluminum Res. Lab., Aluminum Co. of Am., 19i|2. 


I • 


TABLE I 

CHARACTERISTICS OF CYLINDERS 


Cylinder 

(i) 

(in.) 

t 

(in.) 

d i 

(in.) 

d n 

(in.) 

Pi 

(deg) 

PlI 

(deg) 

L 

(in.) 

n 

GALCIT 

65 

10 

0.010 

2.61 

13.05 

15 

75 

i» 

2li 

GALCIT 

25 

15.92 

0.010 

2.53 

22.77 

9 

81 

8 

liO 

PIBAL 

10 

10 

0.012 

3.93 

11.79 

22| 

67 i 

5 

16 


Cylinder 

k» 

€ flat 

€ curved 

€ cr 

(2) 

2w 

(in.) 

Ae ^str 
(sq in.) 

I str r 

(in>) 


G eff 

(psi) 

?i 

^11 


A I 

(lb/in.) 

GALCIT 

65 

5.5 

0.730 x ict 1 * 

2.1*30 x 1CT 1 * 

2.82 x 10’ 1 * 

31.5 x io -li 

0.875 

o.ol*i5 

5.669 x io _lj 

0.2191* x 10"^ 

2.68 x lo 6 

10,180 

255,000 

0.25 

17,500 

GALCIT 

25 

6 

0. 81*7 

1.050 

1.56 

17.0 

1.031 

0.01*29 

5.898 

0.2191* 

3.90 

9,568 

775,ooo 

0.25 

12,331* 

PIBAL 

10 

(3) 

o.5o 

3.0 

3.26 

22.1* 

1.51*7 

0.1596 

22.1*0 

0.8031* 

1.93 

9,862 

88,760 

0.32 

18,190 


•’‘All cylinders are aluminum (E * 10. £ x 10^ psi, G n = 3.9 x 10^ psi, 
v = 0.3, Et from reference 20). u 

2 

In general e must be guessed. In the present investigation e was 
available from experimental results (references 8 and 9). 

not required since € fi a t and € curved are &^- ven in reference 9. 
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TABLE n 

OPERATIONS TABLE FOR THE SIMPLIFIED STRUCTURE 


O 



(1) 

P D2 

(2) 

r B2 

(3) 

^2 

(4) 

n B2 

(5) 

B2 

(6) 

n D2 

(7) 

r C2 

(8) 

C2 

(9) 

n C2 

(10) 

B 1 

(II) 

m t 

1 Cl 

(12) 

r B, 

(13) 

r c, 

(14) 

r 0, 

(1) 

r d 2 

- r >, 

- rr M n 

2a r ? A, 

0 

+ 2a r a t Aj 

0 

0 

rnM, 

-™M n 

+ 2a r a Tl A 1 d 1 

2a r A, 

- 2a r a T A, 

- 2a r a„A,d, 

0 

0 

0 

-2a' r A, 

“^i 

+ 2a r A, 

(2) 

2R 

B2 

0 

-2fr Mj 
- 2Fr Mn 
-4a r ? A t 

0 

2r?i Ml 
- 2rn Mn 
+ 4a r a n A,dj 

2ft M[ 

-2r't Wn 

+ 4a r a, A, 

0 

2a ? r A T 

-2a,a,A, 

- 2a r a n A,d, 

0 

0 

- 2rr Fj 
+ 4a r A, 

-2a 2 r Aj 

0 

(3) 

T D2 

-2a,a r A, 

0 

^7- 

-*S 

-2afA, 

0 

0 

-Os 

-tn»„ 

-2a t a n A,d 1 

-2a,a r A, 

2a, A, 

2a,a„A,d, 

0 

0 

0 

2a,a r Aj 

t'rr, 

-2a,a r A, 

(4) 

2N 

B2 

0 

2nr M , 
-2(Tr Mo 
+ 4a n a r A.d, 

0 

”2iFn Mj 
-2nn u 
-4a A,d, 

-2«„, 

-2n% n 

-4a„a,A,d 1 

0 

-2a„a T A,d, 

2a„a.,A,d, 

2a„A,d' 

0 

0 

2n> t| 

- 4a n a r A.d, 

Sa^Ajd, 

0 

C5) 

2t b 2 

0 

2tr M| 

- 2t? M 
M n 

44a, a r A, 

0 

-2fS M , 

- 2tn M 
-4a,a>,d 

-2fL, 
— 2 n M 

-4a, A, 

0 

-2a,a. r A, 

2a, A, 

2a,a n A l d 1 

0 

0 

2t7 Fl 

-4a,a r A 1 

2a,a r A x 

0 

(6) 

n d 2 

n>M, 

“^M n 

+ 2a n a r A,d, 

0 

-nt M . 

- nt M 

-2a„°a,A,d, 

0 

0 

-nh Mj 

-2w: 

- 2a„a r A,d, 

2a„a,A,d, 

2a ? n A x d, 2 

0 

0 

0 

2a„a r Ajd, 

n7 Fl 

-2a n a r A,dj 

(7) 

2R 

C2 

2af r A, 

2a ? r Aj 

- 2a a, A, 

-2a^„A,d, 

-2a a, A, 

-2 a^Ajd, 

- 2 r~r M 

0 ^ 1 

- 2r :- 0 
-4<A, 

2 n M 
-2 h m ; 

+ 4a r a,A 1 

2rn Mt 

-2fn M 

+4a r a„A 1 d t 

0 

0 

- 2a r Aj 

-2fr F , 
+ 4a' r A, 

- 2a r A, 

(8) 

2 T c 2 

-2a,a r A, 

-2a,a r A, 

2a, A, 

Sa^A.d, 

2a, A, 

2a,a„A,d, 

2fi r«, 
-2tr M , 

+4a,a r A, , 

-2ff^T 
-2»- 
- 4a 2 , A, 

- 2 tn Ml 

“2tn M . 

h -4a^a n A 1 d L 

0 

- 

0 

2ap. r A, 

2t7r, 

- 4a,a r Aj 

2a,a r A, 

(9) 

2 N c 2 

-2a,a r A,d, 

- 2a n a r A t dj 

2a„a, A,d, 

2a n Ajdj 

2a n a,A,d, 

2a„A,d; 

2nr Mi 
2 ht M 

+ 4a n a r A I 4 T 

-2 nt M 
v n 

-4a n a,A l d 1 

-2nn U) 

- 2nn M 
-4<A,V 

0 

0 

2a^. r A,d, 

2nr F , 

-4a n a r A,d, 

2a„a r A,d, 

(10) 

M 

B 1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-fgG-ku) 

”¥D, CkL - s) 

0 

go-') 

0 

(ID 

M 

T CI 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-^(kL-.) 

-f^(j-kLc) 

-go-o 

0 

go-') 

(12) 

R b, 

0 

-2r7 r , 

+ 4^ Aj 

0 

2r7 F , 

-4a r a„A,d, 

2rl F , 

-4a r a,A, 

0 

-2a r A, 

2a r a,A, 

2a r a n A 1 d 1 

0 

-go-') 

-2fr Mj 
-4<A, 
-2Pks/D, 

2a r Aj 

4 Pks/D, 

0 

(13) 

R c. 

- 2a* Aj 

-2<A, 

2a r a T A, 

Sa^V, 

2a A, A , 

2a f a„A I d 1 

~ 2rr r , 

44c£a, 

2r) F , 

-4a r a,A, 

2rn F , 

-4a r a n A,d ] 

go-') 

0 

2a" A j 
4 Pks/D, 

-2r7 
_ 4a r A. 
-2 Pks/D, 

2a r Aj 

4Pks/D, 

(14) . 

V 2 

+ 2a>, 

0 

-2a r a,A, 

0 

0 

r "r. 

-2a r a„A,d, 

-2a r Aj 

2a r a,A, 

2a r a n A,d, 

0 

go-') 

0 

2a ? r A, 
4 Pks/D, 

- rr M , 
-2a>j 

-Pks/D, 
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TABLE IH 

OPERATIONS TABLE FOR SOLUTION WITH ASSUMED DISPLACEMENTS 



0) 

P D2 

(2) 

r B2 

t (3) 

D2 

(4) 

n B2 

t <5> 

B2 

(6) 

n D2 

•P’ 

C 2 

T 

C2 

(9) 

n C2 

(10') 

1 

(1) 

R D2 

" rr M, 

- rr M - 

- 2a r A t 

0 

- rt M 

+ 2?r a i A i 

0 

0 

rn M , 

-rn Mn 

+ a°wM 

2a r A i 

- 2a r a,A, 

- 2a r a„A l d l 

1.0257a* A.- Pr, 

’ 1 F I 

(2) 

2R 

B2 


-2rr„, 

-2rr Uj 

-4<A, 

0 

2rn M , 

-2rn M . 

+ 4a r a n A 1 d, 

2r\,, 

-2rt M „ 

+4a r a l A, 

0 

2 a'A, 

-2^0, A, 

" 2a r a„A,d, 

- 0 64926a* A -006250 Pr f 

(3) 

T Q2 



-‘3m, 

-S 

-2o4 

0 

0 

- tn M 

1 

-2a t a r A, 

2a*A, 

2 a t a rAdi 

-l.0257a t a r A I -»-tr Fi 

(4) 

2N__ 

B2| 




-2nh Ml 

-2nt„, 

-2fit M 

0 

-2a^ r A,d, 

2a -(*l A , d , 

2<Vf 

084926a n a r A,d,+0 06250 nr Fi 

A$) 

2T B2 





-2tt M , 

-2Um d 

-*K a , 

0 

-2a t a r A, 

2 <A, 

^nA.d, 

0 84926a t a r A,+ 0 06250t> F] 

(6) 

N 

D2 






- 

~ nn M, 

-2<A,d, J 

- 2a n<JrM 

2a T,<M, d l 

2 <A,d; 

- 1 0257a n a r Ad + ffr 

F, 

(7) 

2R C2 







r 27F^ 1 

-2r> Mn 

-4a* A, 

2 ^M, 

"2 ft*., 

+ 4a r a,A, 

2r-n Ui 
- 2 rn u . 

+ 4a r a TT A 1 d 1 

- 0. 1 1396a* A.- 0 974 2 6 it 

1 r I 

(8) 

2T C2 








-2ft„, 
- 4a' A, 

— 210m, 

-2tn M , 

- 4<* 1 tt n A,d 1 

Oliaoea^Aj+O 97426 tr F[ 

(9) 

2N C2 









~ 2nh M( 

_2 ""m, 

“ 4a* A,df 

0 1 l396a n a r A,d 1 +0 97426n> Fi 

(10) 

_u 










'(EIL 

-rj- i -f(kL)-094370a r A-l.4765Pr 
L 1 “i 


91592 -S!5ik— 0.44020cm kL—0 47l60kL*inkL + l.356lj 


f— 1 
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TABLE IV 


COMPARISON OF CALCULATED AND EXPERIMENTAL BUCKLING LOADS 


Cylinder 

Method of 
calculation 

Calculated 

p cr 

(lb) 

Experimental 

p cr 

(lb) 

Difference 

(percent) 

GALCIT 25 

Table II 

955 

769 

21*. 2 

GALCIT 65 

Table II 

1670 

1371 

21.8 

GALCIT 65 

Table III 

1730 

1371 

26.2 

PIBAL 10 

Table II 

1*850 

375U 

29.2 


TABLE V 

ASSUMED, CALCULATED, AND EXPERIMENTAL DEFLECTED SHAPES 
[For sign convention and nomenclature see fig. 2] 


Shape 1 r B 

r e 

r D 

m B 

1 

me 


0.06250 

0.56250 

1 

0.22672(1/L) 

0 . 680 l 6 (lA) 


0.031250 

0 . 1*8713 

1 

0.11*170(lA) 

0. 73636 (lA) 


0.015625 

0 . 1*2188 

1 

0.085020(l/L) 

0.76520(lA) 

PIBAL cyl. 10 
(Calculated) 

-0.01691*9 

0.50229 

1 

0. 15086 (1/L) 

0. 80011 (lA) 

GALCIT cyl. 25 
(Calculated) 

0.066323 

0.1*8355 

% 

1 

0 . 093608(1 A) 

0. 82808 (1A) 

GALCIT cyl. 65 
(Calculated) 

-0.16089 

0.3061*0 

1 

-0.21075(lA) 

1.1631(1 A) 

GALCIT cyl. 25 
(Experimental ) 

0.222 

0.667 

1 


— 

GALCIT cyl. 27 
(Experimental ) 

0.0370 

0.201* 

1 

— 


GALCIT cyl. 30 
(Experimental) 

0.0588 

0.1*12 

• 

1 

— 

— 

GALCIT cyl. 35 
(Experimental ) 

0.0893 

0. 1*1*6 

1 
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TABLE VI 

BEAU-COLUMN INFLUENCE COEFFICIENTS 
[^Compression] 


Displacement 

Unit 

Forces 

a) 

f a 

«A 

F B 

m b 

6 A 

Both ends 
fixed 

P ks 

Di 

p 1 — c 
Dl 

P ks 
D 1 

-pifi 

D 1 

m A 


_ p l-c 

D 1 

P s - kLc 
k D-l 

p 1 — c 
D 1 

P Lk - s 
k D-l 

6 b 


p ks 
D 1 

p 1 — c 
D 1 

P ‘ks 

D i 

p 1 — c 
D 1 

m B 


-P 1 - ° 
D 1 

P Lk - s 
k D x 

pi - C 

Dl 

P s - kLc 
k D x 

6 A 

End A 
fixed 

-P — 

tD 2 

-P -~ 

d 2 

P 

tD 2 

0 

m A 

End B 
pin- 
jointed 

-p — 

d 2 

.pi 

d 2 

p — 

d 2 

0 

6 b 

p — 

tD 2 

p — 

d 2 

-p — 

tD 2 

0 

6a 

Both 

ends 

P 

L 

0 

P 

L 

0 

6 b 

pin- 

jointed 

P 

’ L 

0 

P 

L 

0 


^Di = 2(1 - c) - kLs 
D 2 = 1 - (kL/t) 

k 2 = P/EI 


s = sin kL 
c = cos kL 
t = tan kL 



Sign convention 
Forces on constraints 



5,F 


positive as shown 
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TABLE VII 

BEAM-COLUMN INFLUENCE COEFFICIENTS 
[Axial end load P equal to zero] 


Displacement 

Unit 

Forces 

f a- 

m a 

f b 

m b 

6 A 

Both ends 
fixed 

Ml ^ 
Wl M 

C\J 

i — 1 

1 

A EI 

12 — 
L3 

a EI 

m A 

-6 — 

L 2 

L 

A EI 

_2 I! 
L 

6 B 

12 ^ 
L J 

A EI 
6 ? 

-12 51 
L 3 

6 P 

m B 

A El 
■ 6 ? 

_2 EI 
L 

A EI 


6 A 

End A 
fixed 

End B 
pin- 
jointed 


-3 — 

L 2 


0 

m A 




0 

6 B 

3 — 

L 3 


-3 — 

L 3 

0 

h 

Both 

ends 

pin- 

jointed 

0 

0 

0 

0 

6 B 

0 

0 

0 

0 


Sign convention 
Forces on constraints 




I* X 


A 


B 
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TABLE VIII 

BEAM-COLUMN INFLUENCE COEFFICIENTS 
[Tension] 


Displacement 

Unit 

Forces 

(i) 


f a 

m a 

f b 

m b 

6 a 

Both ends 
fixed 

P ks' 
Di- 

(c'_ - 1) 

Dl' 

P ks» 

Dl' 

_p (o' - 1) 

Di 

m A 


r (o' - 1) 
D i 

p/kLc 1 - s'\ 

k \ D i' ) 

r (c - 1) 
Dl' 

p/s* - kh\ 

k l Dp' j 

6 b 


p ks' 

D i 

p (C - 1) 
D l' 

p ks' 
Dl' 

r (c - 1) 
D i 

m B 


p (C - 1) 

Dl' 

p/s* - kL\ 

' k \ D i' j 

p (c - 1) 
D i 

p/kLc * - sA 

k \ D i 1 

6A 

End A 

P k 

-P -i- 

P k 

0 

fixed 

t«D 2 * 

D 2 ' 

t*D 2 * 

m A 

End B 

.pi 

d 2 - 

-P-D- 

d 2 ' 


0 


pin- 

jointed 

p k 

P -1- 

p k 

0 

6 b 

t*D 2 * 

d 2 - 

t * D 2 1 

«A 

Both 

ends 

~ L 

0 

p 

L 

0 

6 B 

pin- 

jointed 

P 

L 

0 

P 

” L 

0 


1 D i » = 2(1 - c 1 ) + kLs 1 
D 2 ' = (kL/t) - 1 

k 2 = P/EI 


s* = sinh kL 
c* = cosh kL 
t * = tanh kL 


Sign convention 
Forces on constraints 


5,F 



P, positive as shovn 
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17S-T4 
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Ring; 

24S-T4 





Cylinders 25 and 65 



Figure 1.- Cylinder characteristics. 
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ym,M 

r,R 


Sign 

convention 





Figure 2 .- Simplified structure for setup of operations tables. 
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Figure 3.- Plot of f(kL) against kL. f(kL) = (-0.91592 

D 1 \ 


sin kL 
kL 


- 0.44020 cos kL 


- 0.47180kL sin kL + 1.3561]. 
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Figure 4.- Determination of buckling load of GALCIT cylinder 25 for two, four, 

six, and eight bays. 



.04 


0 

“nn 

-.04 


-.08 


Figure 5.- Determination of buckling load of GALCIT cylinder 65 by table II. 



Figure 6.- Determination of buckling load of PIBAL cylinder 10 by table II. 
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kL 

Figure 7.- Experimental variation of kL with parameters a and r. 
r/d = 6.32. Value of r given for plotted points. 
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Figure 8.- PIBAL cylinder 82. 
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Figure 9.- Side view of PIBAL cylinder 82 after buckling. 
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Figure 10.- Bottom view of PIBAL cylinder 82 after buckling. 



